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We consider the measuring procedure that in principle allows to avoid the homodyne detection 
for the simultaneous selection of both quadrature components in the far-field. The scheme is based 
■ on the use of the coherent sources of the non-classical light. The possibilities of the procedure are 

^ ' illustrated on the basis of the use of pixellised sources, where the phase-locked sub-Poissonian lasers 

, or the degenerate optical parametric oscillator generating above threshold are chosen as the pixels. 

The theory of the pixellised source of the spatio-temporal squeezed light is elaborated as a part of 
this investigation. 

PACS numbers: 42.50.Dv, 42.50.Lc 



OS 



Oh 1 



Os 
O 



X 



I. INTRODUCTION 



The subject of quantum optics always implies the use of some kind of nonclassical light. Therefore under the 
conditions of the real experiment one of the basic goals is to find a comprehensible source that gives nonclassical light 
with acceptable spatio-temporal properties. For now there are several kinds of sources such as optical parametric 
amplifiers and oscillators generating below threshold (see, e.g., that are non-alternatively used in the practice. 
These devices are preferable mostly because of the nature of parametric interaction itself, which automatically ensures 
a generation of the squeezed light states @, H, 0, [H . 

Together with the obvious advantages of the parametrical sources their applications meet the well-known difficulties. 

■ One of them, for example, is a necessity of the balanced homodyne detection use. In the experimental practice a 
\ preparation of the local oscillator with the acceptable spatio-temporal configuration can be not a simple operation 

■ especially for quantum imaging. So it is reasonable to consider the other alternatives. In this article, the light sources 
are considered to be with quantum fluctuations occurring on high coherent level. We believe that the similar sources 

I . could provide us with new possibilities in comparison with traditional ones. This can be realized by devices such as 
C"*** ' lasers [f| H, Q or by any parametric device generating above threshold [H| • These sources we will call coherent 
sources. As we shall demonstrate further they allow us to avoid an obligative use of the homodyne detection procedure, 
since the coherent component of the radiation itself constitutes for the local oscillator amplitude. 

The development of quantum optics during the last decade has clearly demonstrated that for reasons of any 
applications the greater number of degrees of freedom in the non-classical radiation is a positive and very important 
factor allowing to essentially improve an efficiency of optical measuring procedures. From this standpoint the use of 
the single-mode coherent sources does not seem to be an attractive perspective. However as it was demonstrated in 
Ref. H such conclusion is too hasty because even in this case the radiation outside the high-Q cavity in free space 
turns out to be broadband. 

Having the good temporal statistics the single-mode lasers or the optical parametric oscillation are often ill ac- 
ceptable for the analysis of spatial structures, for example, for the aims of quantum imaging. To adapt them to 
the spatially multi-mode schemes it is suggested to use the so-called pixellised sources that were discussed in Rcfs. 
[13] ■ These sources are formed by putting periodically the pixels (lasers or optical parametric oscillators) on a 
plane surface. That provides us with a possibility to effectively use both the spatial and temporal degrees of freedom 
in distinguish from th e op tical parametric oscillator generating below threshold that was in particular used in the 
dense-coding protocol [12||, where the temporal degrees freedom could not be enable. 

In this article we use the following abbreviations: coherent source (CS); homodyne detection (HD); pixellised source 
(PS); degenerate optical parametric oscillator (DOPO); sub-Poissonian laser (SPL). 
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Figure 1: Illustration of the near- and far-field, / is a focal distance, 2 is a direction of light propagation. 



This article is organized as follows. In Sees. ITll and IIII( two ways (with one and with two sources and correspondingly 
with one and two detectors ) of a direct detection of the quadrature components are considered. In Sees. IIVIVI and 
IVII the measurement procedure is applied for the specific sources, namely for the isolated pointlike phase- locked SPL 
and DOPO, for the PSs on the basis of the SPL and DOPO. 



II. DIRECT DETECTION OF THE QUADRATURE COMPONENTS IN THE FAR-FIELD (SCHEME 

WITH SINGLE COHERENT SOURCE) 

In this article our main aims are coupled with the problems of quantum optics and quantum information that is 
the squeezed and entangled light states are specially important for us. However we would like to emphasize that the 
results obtained in Sees. [TT| and IIIII and devoted to the selection of the quadrature components without the use of 
the homodyne detection technique can have more wide application and are not directly coupled with any quantum 
aspects in the field. 

We shall consider two ways to select the quadrature components presented correspondingly in Fig[T]and FigOl One 
of them is coupled with following the even and odd parts of a photocurrent in the single beam configuration and the 
other with two-beam one. In order to characterize our approach it is important in advance to stress that in distinguish 
from the other approaches we are going to use the so-called coherent sources, which could be formally defined by the 
inequality 

(E N (f)) ^> SE N (r,t), (1) 

where on the left we have an averaged Heisenberg amplitude and on the right its fluctuation. The Heiscnbcrg amplitude 
is written in the form 

E N (r,t) = (E N (?))+5E N (f,t). (2) 

We have a right to say it for both the near and far-fields. Hereinafter the lower index N (F) means the near- (far-) 
field. 

In the limit of a quasi-monochromatic and quasi-plane travelling wave it is convenient to present it in the form [l3j 



E N (z,p,t) = iJp*- e ik ° z ~ S N (p,t). (3) 
V 2e c 

The operator Sn(p, t) is normalized amplitude such that the value (S^ S) takes sense of the mean photon number per 
sec and per cm 2 . Under the propagation in free space the normalized amplitude obeys the canonical commutation 
relations 

S N (p,t),slftfjj\ =5 2 (p-(f)6(t-t'), [s N (i?,t),S N (j?,t')] =0. (4) 

Furthermore, we shall assume that the emission from our sources possess a cylindrical symmetry in the near- and 
far-field: 

(S(p))=(S(-p)). (5) 

By choosing such experimental situation, where a lens is placed at the focal distances / between the source and 
the detection plane the field in front of the detection plane may be considered as a far-field We shall denote the 
corresponding normalized amplitude Sf {p, t) , hence the far- and near- fields turn out to be coupled by the integral 

Mp,t) = -{jJ d 2 P 'S N (f?,t) e-tGf? = ~ § N {$,t), Q=JjP- (6) 



Hereinafter we shall label G(q) as the Fourier image of a function G(p), i.e., 



G(q) = — I d 2 p G{p) e-^P 



(7) 



In particular in Eq. (O we have expressed the Heisenberg far-field amplitude in physical space via the Heisenberg 
near-field amplitude in the Fourier domain with q^Q. 

Let us remind how the wanted quadrature component is selected in the HD technique in the near-field. The 
photocurrent fluctuation operator Si = i — (i^) reads 



6i N (p, t) = p{p) SS N (p, t) + p*5S^ N {p, t) 



(8) 



where (3 is the local oscillator complex amplitude. By choosing the acceptable (3 we can select any quadrature 
component at the detector. For the coherent source it is possible partially to avoid the use of the HD technique 
because, in this case, the mean amplitude (SN(p,t)} takes a role of the local oscillator amplitude. The amplitude 
quadrature turns out to be automatically selected but the phase one remains inaccessible. 

If our measurement procedure is replaced to the far-field, then we shall illustrate that there it is possible to find an 
acceptable approach when the phase quadrature is selected automatically as well as the amplitude one. 

In the far-field the photocurrent fluctuations for the coherent sources are given by 



Si F (p,t) = (S F (p))* 6S F (p,t) + (S F (p)) 6S F (p,t). 



Taking into account ([5]) one can obtain 

5i F {p,t) = 



2- 



(S N (Q))* 6S N (Q,t) + (S N (Q))6S f N (Q,t) 



(9) 



(10) 



Furthermore we shall use the simpler physical situation, where the specific phases are chosen for our sources. In this 
section we require the mean field amplitude in the near field is real, 



(s N (p,t)} = (s N (p,t)y 



(ii) 



In the next section we shall consider the scheme with two sources and for one source we survive the requirement of 
the reality but for the second one we shall foresee the phase shift equal to 7r/2, i.e., our requirement for the second 
source will be (§N(p,t)) = — (Sn(p, £))*■ Then Eq. can be rewritten in the form 



6i F (p,t) 
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(S N (Q)) SS N (Q,t) + h.c. 



(12) 

(13) 
(14) 
(15) 

one can see both quadratures are presented equally in current (|12p . 

Let us vary the direct measurement procedure by investigating not the photocurrent itself but its even and odd parts 
independently. The even an odd parts are proportional to summarized and differential photocurrent respectively: 



Let us introduce the quadrature components Xn(p, t) and Y^(p,t) for the near-field 

S N (p,t) = X N (p,t) + i Y N (p,t), 

and correspondingly in the Fourier domain 

S N (Q,t) = X N (Q,t) + iY N (Q,t). 

Under the assumption 

(Sn(Q)) = (Sn(Q))* = (X N (Q,t)), 



6i±(p,t) = Si F {p,t) ±5i F {-p 1 t). 



Now one can get 



5i+(p,t) 



2- 
Af 
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(X N (Q)) 8X N {Q,t) + h.c. 



5i.(p,t) = 2 ( — ) (X N (Q)) i 5Y N {Q,t) - h.c. 



(16) 

(17) 
(18) 



Figure 2: Scheme of a direct detection of two light entangled beams, BS is a beamsplitter, / is a focal distance. 



One can conclude that the summarized photocurrent is determined by real part of the X-quadrature and the differential 
photocurrent - by imaginary part of the Y-quadrature. On the one hand, it is not the same result that could be obtained 
in the HD technique, where it is possible to select the quadratures as whole. On the other hand, in principle, the 
mathematics provides us with a possibility to renew the full analytical function by knowing only its real or imaginary 
part (the Cauchy-Riemann theorem) . Although this development seems as quite encouraging, nevertheless in the next 
section we want to discuss another measurement procedure, where the quadratures are selected as whole. It is more 
convenient because does not require any additional mathematical treatment. 



III. DIRECT DETECTION OF THE QUADRATURE COMPONENTS IN THE FAR-FIELD (SCHEME 

WITH TWO COHERENT SOURCES) 



Let two independent coherent sources emit two beams (see figH]) presented correspondingly by the slow Heiscnbcrg 
amplitudes Siisr(p,t) and SW(p, £) in the near-field. As usual the corresponding quadratures are introduced by 
relations 



SS mN (Q, t) = SX mN (p, t) + i5Y mN (p, t), 



m = 1,2. 



(19) 



As for our coherent sources we assume that they are perfectly identical but there is a phase shift between them equal 
to 7r/2. In particular this means that the following equalities can take place 



(S 1N (P)) = (Xin(p)), (S 2N (p)) = i(Y 2N (p)), (Y 2N (P)) = (X 1N (p)}, 



and 



(SX 1N (p,t) 5X 1N tf,t')) = (6Y 2N (p,t) 5Y 2N (ff,t')), 
(6Y 1N (p,t) 5Y 1N (f?,t')) = (5X 2N (p,t) 5X 2N (j?,t')). 

After mixing in the near-field on the symmetrical beamsplitter we get two other beams with amplitudes 



Em(p,t) = -y=(S 1N (p,t) + S 2N (p,t)), 



E 2N (p,t) = -^=(S 1N (p,t) - S 2N {p 7 t)). 



(20) 



(21) 



(22) 



These beams can turn out to be entangled if the initial beams were squeezed in the mutually orthogonal quadratures. 
Now we again follow the summarized and differential currents in the far-field that are defined as 



Si±(p,t) = 5i 1F (p,t) ±5i 2F (~p,t), 



(23) 



where i\F and i 2 f are correspondingly the currents on the first and the second pointlike detectors. After non-difficult 
algebraic operations one can obtain 



/ (1 - i) \SX 1N {Q, t) + i 5Y 2N {Q, t)] + h.c \ 



V 1 



5X 2N {Q,t)+iSY 1N (Q,t) 



h.c 

+ h.c 



(24) 



Let us calculate the correlation functions 



(Si±(p,t) 6i ± (f?,t')) = 

= (X 1N (Q))(Xm(Q')) 



5(t-t') 5 2 {Q-Q') + A 



SX 1N (Q,t) 5X 1N {Q',t') :) 
:8Y 1N (Q,t) SY 1N {Q',t') :) 



(25) 



The notation (: ■ ■ • :) means the normally ordered averaging. Under deriving it we have used Eqs. (|21j) staying 
identification of the CS. 

From Eqs. (|25p one can conclude that in the measuring scheme presented in fig. 2 it is possible to select any 
quadrature component by only choosing the summarized or differential current. We shall call this direct measuring 
procedure the plus-minus detection. 



IV. THE PLUS-MINUS DETECTION SCHEME ON THE BASIS OF TWO DOPO OR TWO SPL 

In the previous sections our requirements relative to sources were as follows: they are coherent and cylindrically 
symmetrical. Bellow we consider the specific situations, namely in this section we discuss the plus-minus detection on 
the basis of two phase-locked SPL or two DOPO. 

The quasi-monochromatic radiation with bounded aperture inside the cavity at the output mirror can be formally 
described by the Heisenberg amplitude 



E s (z = L,p,t) =i\j^ e~ lUJot f(p) a(t). (26) 

Here p is the transverse spatial vector, 2 is a longitudinal co-ordinate, L is the cavity perimeter, luq is the mode 
frequency. The normalized amplitudes a(t) and a^(t) obey the canonical commutation relation 

[a(t).a + W] =1 (27) 

and (tf(t)a{t)} gives the mean photon number inside the cavity. 

The transverse size of the mode is formally bounded by the function f(p) that is normalized such that 

d 2 p \f(p)\ 2 = 1, d 2 p = dxdy. (28) 

For the sake of simplicity we shall discuss here the Gaussian mode, then 



1 . ( P 



■2 



The value w$ takes a sense of the light spot size. 

The Heisenberg amplitude d(t) describes the intracavity field. Certainly, its time evolution depends on the specific 
physical processes inside the cavity. For the coherent field 

a(t) = (a) + Sa(t), (a) > 5a. (30) 

The near-field at the output of the cavity with amplitude E^{p, i) is formed as a linear combination of the field leaving 
the cavity for a propagation in free space and the multi-mode vacuum field reflected from the output mirror outside 
the cavity 

E N (p,t) = VTE s {z = L,p,t) -VRE vac (z = -0,p,t). (31) 

Here R and T are correspondingly the reflection and transmission coefficients. We put there are no losses of the field 
under passing the mirror, i.e., R+T=l. 

Taking into account ([3]) and (|26|) . we can rewrite (|3ip via the normalized amplitudes in the form 

S N (p,t) = yfcf(p) a(t) - S vac (p,t), (32) 



where 



cT 



and R « 1 , 



(33) 



Here the normalized amplitude S vac (p,t) as well as £>N(p,t) obeys free-space commutation relations (jl]). In order to 
ensure this we should put 



^[a(t),a\t')} = \a(t),Sl c (t')] + \s vac (t)A f (t') 



(34) 



Because we want to consider the ztdetcction we need to foresee two sources and correspondingly we have to discuss 
the amplitudes S riVac (p, t), S r ^(p, t), and a r (t), where index r = 1, 2 indicates the source. We shall assume that both 
sources are perfectly identical with each other, and the parameters k and f(p) are the same for both. However it is 
suggested that there is the phase shift between them such that 



{a 2 ) = i\Jn, 



and equalities (|2Tj) take place. 

Applying these considerations to our sources, it is not difficult to get the quadrature components 

8X r ,N(p,t) = \/k f(p)Sx r (t) - X r>vac (p,t), 
8Y r>N (p,t) = Vk f{p)6y r {t) -Y r>vac (p,t), 
(X hN (p)) = (Y 2 , N (P))=V^f(p), 
(X 2 ,n(p)) = (Y 1 , N (p))=0, 

where Sx r (t) and Sy r (t) are the quadrature components of the intracavity field for both sources 

8a r (t) = 8x r (t) + iSy r (t). r = 1, 2. 

Passing on to the Fourier domain, one can obtain 

Sx it N(Q,t) 
5yi,N{Q,t) 



/k f$8xi(t) - Xi tVac (Q,t), 

^ fn *2?i(*) - yi,vac(Q,t), 



Here 



and for the Gaussian mode 



(xi,n(Q)) = (V2,n(Q)) = V^ri fn- 



f$=ll d* P f(p)e-*QP 



(35) 

(36) 
(37) 
(38) 
(39) 

(40) 

(41) 
(42) 
(43) 

(44) 
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(45) 



Now we have a complete information about our sources to calculate a final correlation function. Substituting (|41j) -(|4"3" ]) 
to Eq. |2"5|) . one can obtain 



(Si±(p,t) 8i±{ff,i?)) = Ann 



f of & 



8(t-t') S 2 (Q-Q') + 4k 



(: fcci(t)azi(f) :) 
(: 8m(t)mt') ■■) 

Let us consider the spectrum of the spatio-temporal correlation function determined as 



2n 

A? 



(46) 



- lim — 



T/2 pT/2 



T/2 J -T/2 



dtdt' [ [ d 2 pd 2 p'(5i ± (p,t)5i ± (f/,t')) e itt(t-t') e iq(p-f?) (47 ) 



Substituting (|46[) to (|4T[) after all integrations one can obtain the spectra in the explicit form. For DOPO it can be 
read via the physical parameters: 



Ann 



1 ± e ~wW 



- A/ 1 

™0 = 7y ■ 



(48) 



Under calculations of these formulas we have used the expressions for the correlation functions (: 5xi(t)5xi(t') :) and 
(: Syi(t)Syi(t') :) obtained in Ref. fill ]. Here pith > 1 is a factor determining a relationship between the pump power 
and the threshold one. 

By the same way we could find the spectra for phase-locked SPLs in the scheme of plus-minus detection. Using the 
results obtained in Ref. [1, 0] one can read 



Ann 



( 



-2 2 

-w q 



n 2 



2k 2 

V n 2 /4 + n 2 J 



(49) 



It is not difficult to check that the measurements in the HD technique in the near-field will give exactly the same 
results. The last formulas exhibit the possibility of the plus-minus detection in the cases of the simplest single-mode 
sources. One can see that it is possible to follow both the squeezed and stretched quadratures simultaneously without 
the use of the HD technique. 



V. PIXELLISED SOURCE: ARRAY OF THE COHERENT POINTLIKE SOURCES 

Let us consider more complicated sources that we call the PS. Let N 2 pointlikc identical sources be placed period- 
ically on the plane such that the linear distance between the adjoining pixels is much more than the beam spot each 
of the individual pixels 

I > w . (50) 

This requirement allows us to say that all sources are perfectly independent from each other. Simultaneously we 
assume that all the lasers or DOPO used as the pixels are perfectly synchronized by means of the pump field for the 
DOPO or by means of the special synchronizing external field for the lasers [f| 0]- 

In order to construct theoretically the full beam from the PS we can use the linear superposition principle and 
derive the full Heisenberg amplitude as a sum of the individual amplitudes from each of the pixels, 



E{r,t) = Y J E m {r,t). (51) 



The positions of the pixels on the surface is determined by the vectors 

p m = lm, m x ,m y = 0,±l, ±2,..., ±(N-l)/2. (52) 



In Eq. (|32|) we have written the explicit form for the slow Heisenberg amplitude in the case of the single pixel. It is 
not difficult to generalize it on the case of the PS 

S(p,t) = f(p- p m ) a m (t) - S va ,c(p,t). (53) 

Pm 

The function f{p — p m ) describes the m-beam aperture. Since we choose I 3> wq, then all the pixels are independent 
from each other and we have a right to require that 

\f(p-p rn )\ 2 cfp=l (54) 



and 



[a m (t),al{t)] = 8„ 



(55) 



The operators a m are the normalized amplitudes describing a single mode field inside each of the individual cavities. 
As before we consider the coherent fields that is 



a m (t) = (a m ) + 8a m (t), (a m ) > 8a m {t). 



(56) 



Because under plus-minus detection we have to foresee two independent PSs we are under obligation to specify 
belonging the m-th pixel to one of these PSs that is to add additional index: a m — > a r , m , where r = 1,2. 
Now it is not difficult to generalize formulas (jlTj) - (|4"3"f 



8X r>N (Q,t) = ^/qXI^^vW - X nvac (Q,t), 
5Y riN ($, t) = y/Hf^ ^ m Sy r , m {t) - f^ ac (Q, t), 

Here 8x r ^ m (t) and 8y r ^ m (t) are the quadrature components for the intracavity single- mode field 

5a T}m = 8x rtTn {t) + io~y r}m {t) 

and 



(57) 
(58) 
(59) 

(60) 
(61) 



Evaluating the series in Eqs. (|6ip one can get it in the explicit form, 

sinQJN/2 sinQ y lN/2 
smQ x l/2 sinQyl/2 ' 



(62) 



where Q x and Q y are the components of vector Q. On the basis of Eqs. (|5T|) - ((5"9"|) the correlation functions can be 
read 



(Si±(p,t) St ± (f?,t')) = 



1^) f Q A Q 



8{t-t')5 2 {Q-Q')+An 



[: 8x r (t) Sx r (t') :) 
(: Sy r (t) 8y r (t') :) 



(63) 



Passing to the Fourier domain, one can obtain again for the DOPO and SPL Eqs. (|48[) - fl9")) . where we should do the 
following exchange 



-2 -2 
e -w q 



1 



E 



w (q + Q m - Q n ) 



2tt 



(64) 



In Fig. [3] one can see that the photocurrent spectrum exhibits the periodical structure in dependence on spatial 
frequency q that reflects directly the periodical geometry of the PS. The reduction of the shot noise is more effective 
for bigger N, where the bigger number of holes with nearly perfect noise reduction appears. 



VI. DIRECT DETECTION IN THE SINGLE BEAM 



Let us return to the measuring procedure presented in FigfTJ where we follow the even and odd parts of current 
under the detection in the beam from the single PS. 
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Figure 3: Spatial spectrum of the summarized (for SPL) and differential (for DOPO) photocurrent variances; q y = 0, f2 
0, l/wo = 10 in arbitrary units, N = 7 (a), N = 99 (b). 



Taking into account the formulas from the previous section it is possible to get the explicit expression for the pair 
correlation function of the summarized and differential currents in the form 



+4/v 



5x r (t) 6x r (t') :) 
6y r (t) 6y r (t') :) 



fd fd< \ ^d+d' ^ ^-6-3' 



Using again the definition |17|) it is possible to get the spectra for the DOPO and SPL in the form Eqs. 
(HU), where we should do the following exchange 



-2-2 

-w q 



2N 2 



E 

Q m Qr, 



-w 2 (q+ \{Q m - Qn)j (q + Q m - Q., " 

e V / ± e V 



(65) 
and 

(66) 



The formula exhibits a set of resonant Gaussian picks or holes on the level of the shot noise (see figs f4] and [5]). Each 
of the picks has a width 



1 2tt 
= — = — w , 
w A./ 

and is centered on the frequency Q m or Q m /2 that matches the pixel positions: 

n - 2w z 

tym — ^ n Prn ■ 



(67) 



(68) 



The distance between the adjoining picks along the x- or y-axes d and the full actual frequency range D are given by 

(69) 



2ir 2tt 
d= — I and D = — IN. 
A/ A/ 



Let us remind that values wq, I (wq <C I), and N, respectively, corresponds to the linear size of the pixel, a distance 
between the adjoining pixels, and the pixel number along one of the transverse directions. 

We want to follow the reduction of the shot noise, i.e., the summarized current in the case of the SPL pixels and 
the differential one in the case of the DOPO pixels. In (figH|), an interchange for the holes on the main frequencies 
(in case l/wo = 10 they are frequencies with numbers divisible by 10 and zero frequency) and half less deep holes on 
the additional frequencies (here with the numbers divisible by 5) occurs. 

Both sequences of the holes decline to the shot noise level on the frequencies far from zero. The effective hole 
number is determined by the parameter D (Eq. ([fj9^1 ). One can see that for small N (see fig'Hh); complete reduction 
of the shot noise takes place only around zero frequency q x . For the greater number of N (see fig©) , the interval 
which damping of the holes in the shot noise level occurs on gets larger. Therefore the well reduced fluctuations can 
be observed at the wider range of q x . 

The situation when the DOPO is chosen as the pixels (figEJ) is unfavorable since deeps at the shot noise level 
dependency are alternating with picks. Moreover, the depth of the holes is limited by half of shot noise level. 
Therefore the perfect suppression of the shot noise at this case is simply impossible. 
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Figure 4: Spatial summarized photocurrent variances spectrum in the far field, in case when SPL is chosen as the pixel, 
q y = 0, f2 = 0, l/w = 10 in arbitrary units, N = 7 (a), N = 99 (b). 
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Figure 5: Spatial differential photocurrent variances spectrum in the far field, in case when DOPO is chosen as the pixel, 
q y — 0, SI = 0, 1/wq = 10 in arbitrary units, N = 7. 



VII. CONCLUSION 



In this article wc purpose two main aims. First of all we wanted to find the measuring procedure that could be in 
some cases an alternative to the balanced homodyne detection. The latter is an universal approach, which allows to 
select any quadrature component of the field by mixing on the detector with the field from the local oscillator. For 
that we have to choose the acceptable phase and the acceptable spatio-temporal configuration of the local oscillator 
amplitude. Often especially in quantum imaging, where the spatial configuration can be extremely complicated, this 
task turns out to be not very simple. Then it would be nice to find an approach, where these obstacles could be 
graded to some way. 

To our mind we have discussed the quite adequate scheme that allows to follow both squeezed and stretched 
quadrature components simultaneously without the use of the HD technique. For that we need first to use the 
so-called coherent sources and second to follow not current itself but the summarized and differential combinations 
currents under the detection of two entangled beams in the far-field. 

The second goal was as follows. We constructed the theory of so-called pixellised source of the spatio-temporal 
squeezed light. Wc assumed to use the phase-locked sub-Poissonian lasers and the optical parametric oscillators 
generating above threshold as pixels. On this basis wc have examined our measuring scheme. 
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